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ABSTRACT 

In addition to producing a strong gravitational signal, a short gamma-ray 
burst (GRB), and a compact remnant, neutron star mergers eject significant masses 
at significant kinetic energies. This mass ejection takes place via dynamical mass ejec¬ 
tion and a GRB jet but other processes have also been suggested: a shock-breakout 
material, a cocoon resulting from the interaction of the jet with other ejecta, and vis¬ 
cous and neutrino driven winds from the central remnant or the accretion disk. The 
different components of the ejected masses include up to a few percent of a solar mass, 
some of which is ejected at relativistic velocities. The interaction of these ejecta with 
the surrounding interstellar medium will produce a long lasting radio flare, in a similar 
way to GRB afterglows or to radio supernovae. The relative strength of the different 
signals depends strongly on the viewing angle. An observer along the jet axis or close 
to it will detect a strong signal at a few dozen days from the radio afterglow (or the 
orphan radio afterglow) produced by the highly relativistic GRB jet. For a generic 
observer at larger viewing angles, the dynamical ejecta, whose contribution peaks a 
year or so after the event, will generally dominate. Depending on the observed fre¬ 
quency and the external density, other components may also give rise to a significant 
contribution. We also compare these estimates with the radio signature of the short 
GRB 130603B. The radio flare from the dynamical ejecta might be detectable with 
the EVLA and the LOFAR for the higher range of external densities n > 0.5cm -3 . 

Key words: gravitational waves—binaries:close—starsmeutron—gamma-ray 

burst :general 


1 INTRODUCTION 


A binary neutron star (ns 2 ) merger is one of the most 
promising targets of ground-based gravitational-wave (GW) 
interferometers, such as Advanced LI GO, Advanced Virgo , 
and KAGRA dThe LIGO Scientific Collaboration! 120141 : 
lAccadia et al.l 1201 ll : lAso et alJ 20131 1. The expected even t 


rate of ns 2 mergers is 0.4 - 400 yr -1 jAbadie et alj|20ldl . 
Most of these events will be just above or just below the de¬ 
tection threshold. Observations of an electromagnetic coun¬ 
terpart will confirm the validity of these GW signals, increas¬ 
ing signific antly the potential detect ion rate and detection 
confidence dKochanek fc Piranlll99 j f. In addition, an elec¬ 
tromagnetic counterpart will enable the localization of the 
sources and the identification of their host galaxies and their 
redshifts, enhancing significantly the potential information 
from this event. Finally any electromagnetic counterpart will 
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provide invaluable information on the physics of the merger 
process. 

The detection horizon distance will extend up to a few 
hundred Mpc. The size of the GW-sky localization error 
box will depend on the number of detectors used and be¬ 
tween a few tens and thousands Milky- Way size galaxies wil l 
reside within t his error box (see e.g., iNissanke et all 1 201 ll ; 
iFairhurstl 1201 il l . Follow-up observations will be a challeng¬ 
ing task, even for a search limited to these galaxies. Clearly, 
a good understanding of the expected electromagnetic sig- 
nals is essential to detect an electromagnet i c counterpart 


(lMgtz ger_fc_^erge3|20lj; Nissanke et al.ll2013l:lKanner et ahl 


120131 : Kasliwal fe Nissankd 20141 ; Bartos et al.l 2014) 1. 

Ns 2 mergers have been recognized as the possible pro¬ 
genitors of short g amma-ray bursts (GRBs) (lEichler et al.l 
119891 : iNakarl 120071 1 and short GRBs and their afterglows 
are one of the most attractive electromagnetic counterparts 
of GW events. However, GRBs and their early afterglows 
are believ ed to be highly beamed with a half-opening angle 


9j « 10° jFong et al.l 2014) 1. This results in about 5% prob- 
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Table 1 . Summary of mass ejection in different forms. 



Mass [Mg] 

Kinetic Energy E [erg] 

Fiducial E [erg] 

Average /3T 

Fiducial /3F 

Average Y e 

Reference 

Dynamical ejecta a 

10~ 4 - 10~ 2 

10 4a - 10 &i 

5 ■ 10 8u 

0.1 - 0.3 

0.2 

0.01 0.4 

[i] 

GRB jet 

< 10-8 

io 47 - to 50 - 5 

10 48 

> 30 

- 

- 

[2] 

Cocoon 

10~ 6 - 10" 4 

10 47 - to 50 - 5 

10 48 

0.2 - 10 

0.3 

- 

[3] 

Shock breakout 

0 

1 

o> 

1 

0 

1 

IQ 47 _ 10 49.5 

to 48 - 5 

1 

1 

- 

[4] 

Wind^ 

10- 4 - 5 ■ 10- 2 

10 47 - 10 5 ° 

10 5 ° 

0.03 - 0.1 

0.07 

0.2 - 0.4 

[5] 


References; 

[1] iGorielv et al ] d201lh ; iKorobkin et al.l d 201 2b ; I Hotokezaka et al.l d2013h ; iBauswein et al.l d2013h ; iRosswoel d2013h ; I Piran et al.l d2013h ; IWanaio et all d2014l) , 

[2] | N akarl j2007|) ; IWanderman 8z Piranl d2014|) , 

[3] iNagakura et al.l d 201 dl ') ; iMurguia-Berthier et al.l d 20141) , 

[4] iKyutoku et all d2014| ) ; I Metzger et al.| d2015|) , 

[5] iDessart et al.l d2009l) ; iFernandez fc Metzger] d2013l) ; iMetzger &c Fe rnandez ( 2 0141) ; I Per ego et al.l d2014f) ; I Just et al.l d2014l) ; I Fernandez et all d2015l) . 

a This component is composed by the tidal tail and shocked component. The main difference between them is the value of Y e . The tidal 
component can have a lower Y e . 

^ This component includes two cases depending on whether the remnant is a black hole or a neutron star. In the former the wind is just 
from the surrounding disk while in the latter it arises from the neutron star as well. The wind from the remnant neutron star has a higher 
value of /5T and Y e . Note that the fiducial value used is an optimistic one. 


abilit MJ to detect a s hort GRB in coincid e nce with t he GW 
signal llSchutzlIioill : iNissanke et ahll201ll : ISetoll2015h . Even 
if the viewing angle is larger than the jet opening angle, off- 
axis afterglows, called orphan afterglows, can be observed at 
late times when the relativistic jet slows down and its emis¬ 
sion is less beamed. As the frequency of the peak flux of 
afterglows decreases with time, at sufficiently small view¬ 
ing angles the off-axis afterglows in the optical to radio 


bands can be a good potent ial candidate of electromagnetic 
counterparts to GW eve nts (Ivan Eerten fc MacFadvenlhoill : 
IMetzger fe Bergeill2012fl . 

In addition to the GRB beamed emission and its late 
more isotropic orphan afterglow, electromagnetic waves will 
be emitted quasi isotropically at different stages from ma¬ 
terial that is ejected during the merger. Most notable 
one is a macronova (also called kilonova), an optical- 
infrared transient driven by the radioactive decay of the 


Metzger et al. 

120101: Kasen et al. 

20131: Barnes & Kasenl 

20131.: Tanaka & Hotokezaka 201S : 

Grossman et al. 2014f). 


Recently, the Hubble Space Telescope detected a near 
infrared bump at 9 days after the Swift s hort GRB 
130603B dTanvir et al] I20I3I ; iBerger et al.l EoTiil i. which is 
consistent with the theoretical expectation of macronovae. 
While the identification hinges on a single data point if cor¬ 
rect this is the first observational evidence for a significant 
mass ejection with a high velocity from a ns 2 merger. 

Synchrotron radiation of electrons accelerated in shocks 
formed between the (mildly) relativistic ejecta and the in- 
terstclla r medium (ISM) is a second electromagnetic coun¬ 
terpart dNakar fe Piranlboilh . This emission can last up to 
a few years and peaks in the radio band. All the ejected 
material will contribute to this emission, but different com¬ 
ponents with different velocities will contribute at different 
timescales, at different frequencies, and at different intensi¬ 
ties. The rise time and the peak flux depend on the density 


1 Note that these estimates take into account that the GW hori¬ 
zon is larger in the direction of the GRB hence the chance of 
coincidence with a short GRB is larger than the beaming frac¬ 
tion. 


of the ISM surrounding the merger but for modest densi¬ 
ties the radio signals can be observed up to the detector 
horizon (iPiran et al.ll2013l f . 

Recent studies have shown that mass ejection from 
mergers is driven by several different processes. The most 
robust one, that appears in numerical merger simulations, 
is the dynamical mass ejection. If mergers are accompa¬ 
nied by GRBs then clearly relativistic GRB jets are an¬ 
other component. Other mass ejection mechanisms that 
have been proposed are merger shock-breakout material, vis¬ 
cous/neutrino/magnetically driven winds, and a possible co¬ 
coon that forms when the jet propagates within the other 
components of the ejecta. 


The different components have different masses and 
kinetic energies. Their characteristics depend on the na¬ 
ture of the progenitors, in particular on their rela¬ 
tive sizes, on the nature of the merger remnant, which 
could be either a black hole or a massive neutron 
star (MNS), as well on the, unk nown yet, neutron star 
matter equation of st ate (see, e.g., iHotokezaka et ahll2013l : 
IBauswein et alj 120131 ) . The different components will in¬ 
teract with each other and these interactions will affect 


20141; Murguia-Berthier et al.1 20141), nossiblv producing 

electromagnetic signatures (e.g., 

Bucciantini et al. 

2012; 

Zhang] i2013l: IMetzger & Pirol .2014 

: Nakamura et alj 

2014; 

Rezzolla & Kumarl 20141; Ciolfi & Siegel 20141: Kisaka et al. 

2014)). 

In this paper we examine the long lasting radio 

emis- 


sion arising from the different components of the ejecta. The 
structure of the paper is as follows: We summarize in Sec.O 
the properties of the different components of the ejecta. In 
Sec.[3]we calculate the expected long-lasting radio flares pro¬ 
duced by the interaction of the different components of the 
ejecta with the ISM. In Sec. [4] we compare these estimates 
with the radio signature of the short GRB 130603B. Finally, 
in Sec. 0 we summarize our results and their possible im¬ 
plications on the detection of radio signals accompanying 
mergers. 
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Mass ejection and radio signals from NS mergers 3 



Figure 1. Left panel: the kinetic energy and the four-velocity of the different components of the ejecta. Also marked are the deceleration 
timescales of Eq. m assuming an external density of 1 cm 3 . The star in each component shows the fiducial model. Right panel: a 
schematic picture of the morphology of th e different components o f the ejecta on the meridional plane. The distribution of the dynamical 
ejecta is taken from a merger simulation IIHotokezaka et al.ll20 1 .Ti . Other components are added schematically. 


2 DIFFERENT COMPONENT OF EJECTA 

AND THEIR PROPERTIES 

As material is ejected in different processes the different 
components will have different masses, kinetic energies, ve¬ 
locities, and electron fractions. The first three quantities 
determine the radio flare signals while all four are impor¬ 
tant for macronova estimate. Table [T] summarizes the values 
of these quantities as taken from the recent literature. The 
properties of the different components of the ejecta are also 
shown in Fig. |T] The left panel of the figure depicts the pos¬ 
sible range of the kinetic energy, E, and the four velocity, 
T/3. Here T is a Lorentz factor and f5 is a velocity in units 
of the speed of light c. Also shown in the figure are the 
deceleration timescales due to the interaction with the ISM, 
which are discussed later. This timescale gives the character¬ 
istic peak time of the radio flares from each component. The 
right panel of the figure shows schematically the expected 
morphology of the ejecta. 

In the following, we briefly describe the properties of 
the different components. In each case we focus on the to¬ 
tal mass, energy, and the corresponding velocities. We also 
mention the expected distribution of energy as a function 
of velocity, which is essential in order to estimate the ra¬ 
dio flares from these components. For completeness we also 
mention the electron fraction Y e . This is not needed for the 
radio estimate but it is a critical quantity that determines 
the composition of the ejected material as well as the heating 
rate that is essential for macronova estimates. 


2.1 The dynamical ejecta 


Gravitational and hydrodynamical interactions produce the 
dynamical ejecta. In many senses it is the easiest to calcu¬ 
late and as such it is the most robust ele ment. It was investi¬ 
gated using Newtonian simulations (e.g . Davieset_alJ 1994 
iRuffert et al.lll997l ; iRosswog et al.lfl 9991 : Rosswoe 2012 ) and 
using general relativistic sim u lations (e.g., Oechslin et al.l 
120011 ' IHotokezaka et all 12014 IBauswein et al. 20l4 ). Ac¬ 


cording to these numerical simulations, the mass and ki¬ 
netic energy of the dynamical ejecta are expected to be in 
the range 10 -4 < M e j < 10 _2 Mq and 10 49 < E < 10 51 erg, 
respectively. The median value of E in the general relativis¬ 
tic simulations is a few times 10 5 ° erg. The properties of the 
dynamical ejecta are as follows. 

The tidal ejecta. A fraction of the material obtains suffi¬ 
cient angular momentum and is ejected via tidal interaction 
due to non-axisymmetry of the gravitational forces. This 
matter is ejected even before the two stars collide with each 
other and it lasts as long as the gravitational field is not 
axisymmetric (about 10 ms after the merger in the case 
that the remnant is a MNS). This tidal component is mostly 
ejected into the equatorial plane of the binary within an a n- 
gle about 20° (see e.g., Fig. 17 in IHotokezaka et akll20l4 l. 

The electron fraction of the dynamical ejecta and 
the resulting nucleosynthesis have been studied in the 
literature (e.g. iGorielv et al.l 1201 ll : iKorobkin et al] 120121 ; 
IWanaio et al.1 120l4 ~l. The tidally ejected material has ini¬ 
tially a low electron fraction Y e -C 0.1 as this matter 
does not suffer from sho ck heating and neutrino irradia¬ 
tion (IWanaio et al.l l2014h . This is particularly important 
concerning the possibility that this is the source of heavy 
(high atomic number) r-process nuclides, but it is not so rel¬ 
evant for our discussion that is concerned mostly with the 
radio flare. This fraction can increase by electron neutrino 
absorption or by positron absorption. The tidal component 
ejected at late times has higher Y e values. 

The shocked component. A shock is formed at the in¬ 
terface of the merging neutron stars. The shock sweeps up 
the material in the envelope of the merging neutron stars. 
Furthermore, a shock is continuously produced around the 
envelope of a remnant MNS as long as the MNS has ra¬ 
dial oscillation. As a result, a fraction of the shocked ma¬ 
terial obtains sufficient energy and is ejected from the sys¬ 
tem. Recent general relativistic simulations show that this 
component can dominate over the tidal component in the 
case o f a nearly equal mass binary (e.g., IHotokezaka et ahl 
12014 IBauswein et al] I20l4 f. The shocked component is 
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4 K. Hotokezaka and T. Piran 


ejected even in the direction of the rotation axis of the 
binary. The average electron fraction of the shocked com- 
ponent s is relatively large compared with that of the tidal 
ejecta dWanaio et alj|2014f l. It may be as large as Y e ~ 0.2 - 
0.4 and it will result in a different nucleosynthesis signature. 

We take the velocity distribution of the dynamical 
eje cta from the resu lt of a numerical relativity simulation 
of iHotokezaka et ahl (l2013h for a 1.4- IAMq ns 2 merger 
for the case of APR4 equation of state. The energy dis¬ 
tribution of this model can be approximately described as 
E( ^ p) oc /3~ 0,5 with a cut off at p ~ 0.4 and an average 
velocity is P ~ 0.2, where E(J^ p) is the kinetic energy with 
a velocity larger than /?. Note that it is not clear whether the 
cut off at P ~ 0.4 is physical or that it arises just because it 
is difficult to resolve such a small amount of fast material in 
the numerical simulations. For our fiducial model, we use a 
total kinetic energy of 5 x 10 50 erg. 

The relativistic shock-breakout component. When the 
shock breaks out from the neutron star surface to the ISM, 
it is accelerated and a fraction of the shoc ked component 
can h ave a relativistic velocity with PT > l. lKvutoku et al .1 
d2014lf showed analytically that the kinetic energy of the 
relativi stic ejecta can be_~ 10 4 ' 10 49,5 erg. More re¬ 

cently, I Metzger et al.l (|2015i l found that there is a mildly 
re lativistic component w ith P > 0.8 in a merger simulation 
of lBauswein et aD d2013l l. This fast component is likely re¬ 
sulted from the acceleration of a shock emerging from the 
neutron star surface. They found that the mass and kinetic 
energy of the fast component with P > 0.8 are ~ 10 _5 Mq 
and ~5 x 10 48 erg, respectively. Because of this large veloc¬ 
ity, the radio signature of this component would be differ¬ 
ent from the slower material. We denote this component as 
a “shock-breakout material” and we consider it separately 
from the sub-relativistic dynamical ejecta. 

Here, we assume that the kinetic energy distribution 
of this component is a simple power-law, E (^ p r) = 
10 48 ’ 5 ( pT)~ a erg as the fiducial model. The value of a varies 
from 1.1 f or ,'ffi 1 to 5.2 for /3F <C 1 dKvutoku et al.lf2014l : 
iTan et ahl 1200il l. We set a to be 3, which is valid around 
pr ~ 1 and we take into account only the fast component 
with pr JS 1. 


2.2 The ultra-relativistic jet 

If ns 2 mergers are progenitors of short GRBs, they involve 
relativistic jets. Assuming the kinetic energy of the jet as 
the gamma-ray energy of the prompt emission, the energy of 
the jet can be estimated from the observed GRBs. The mini¬ 
mal and maximal values of the observed isotropic-equivalent 
gamma-ray energy for non-Collapsar short GR Bs are 2 x 
10 49 erg and 4 x 10 52 erg, respectively dWanderman fe Piranl 
l2014l l . Taking into account the averag e value of the mea sured 
jet-half opening angles of 0j « 10° dFong et aid 1201 il l, the 
kinetic energy of a relativistic jet is in the range of 10 4 ' - 
3 x 10 5 ° erg. The luminosity function is rather steep and 
there are more weak GRBs than strong ones. Hence we con¬ 
sider here a fiducial GRB jet with a kinetic energy of 10 48 erg 
and a jet-half opening angle of 10°. 


2.3 The wind from the merger remnant 


Some of the debris of the neutron stars form an accre¬ 
tion disk that surrounds the central remnant. The mass 
of this accretion disk is estima ted to b e in the range of 
10~ 3 < Mdisk 5-, 0.3 (see e.g., S hibata fe Taniguchil2006l : 
IRezzolla et alj|2010l : IHotokezaka et al.ll2013l l . This accretion 
disk produces an outflow driven by viscous and neutrino 
heating. The properties of this outflow depend on the cen¬ 
tral object as follows. 

A black hole with an accretion disk. The wind from an 
accretion disk surrounding a black hole has been explored , 
in the context of mergers, by Fernandez fe Metzger! (I2013lj : 
ljust et ahl iH3) ; iFernandez et alj ( 2015T ). The disk is ex¬ 
tremely dense and the accretion rate is huge. Initially it is 
opaque even for neutrinos. After ~ 0.1 - 1 s from the onset 
of the merger, the density and temperature of the accretion 
disk decrease and neutrino-cooling becomes inefficient. As a 
result, a fraction of the material is ejected isotropically due 
to the viscous heating in the accretion disk. The amount of 
ejected material is about 5 - 20% of the initial disk mass 
depending on the a-viscosity parameter and on the spin pa¬ 
rameter of the black hole. Increasing these parameters, the 
fraction of the ejected mass to the initial disk mass increases. 
The average velocity of the ejecta is 0.03 < P < 0.05 and 
the expected kinetic energy of the outflow is in the range 
of 10 47 < E < 10 5 ° erg. The average electron fraction is 
Ye ~ 0.2 - 0.3. 

A neutron star with an accretion disk. The wind from a 
neutron star with an accretion disk can be divided into three 
parts; a neutrino-driven wind from the remnant neutron star 
itself, a neutrino driven wind from the accreti on disk, and a 
viscous driven wind from the acc r etion disk dD essart e t al.l 
120091 : iMetzger fe Fernandezl 120141 ; IPerego et alJ 2014 ). The 
neutrino-driven wind expands into relatively high latitudes 
and it has a larger velocity P < 0.1 and a higher elec¬ 
tron fraction Y e ~ 0.4 than those of the viscous-driven 
wind. The amount of material ejected by the neutrino-driven 
wind depend s on t he lifetime of the central neutron star. 
IPerego et alj J2014l l showed that the ejected mass is more 
than 3 x 10 _3 Mq with a velocity P ~ 0.06 - 0.9 at 100 ms 
after the merger in the case o f an initial disk mass of 0. 17Mq . 
IMetzger fe Fernandez! (|2014l l also showed that the amount 
of the ejected mass and average velocity are M e j ~ 10~ 3 Mq 
and P ~ 0.05 at 100 ms with an initial disk mass of 0.03Mq. 
When a MNS does not collapse into a black hole, about 20% 
of the initial disk mas s may be ejected as neutrino driven or 
viscous driven winds (IMetzger fe Fernandedl2014l l. 

In the following, we take the kinetic energy of 10 5 ° erg 
with a single velocity 0.07c for the wind from the merger 
remnant as the fiducial model. Note that these values cor¬ 
respond to the most optimistic case. As shown later even in 
this case, the expected radio signals are very weak. 


2.4 The cocoon 

The interaction of a GRB jet with the pre-ejected mate¬ 
rial such as the dynamical ejecta or the wind along the ro¬ 
tation _axis_wouId_OToduce a hot cocoon surrounding the 
iet (iN agakura et al] 1 20141 : iMurguia-Berthier et al.l 120141 : 
IRezzolla fe Kumarll2014 1. After the jet emerges from the ex¬ 
panding ejecta, the cocoon will break out from the surface 
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of the ejecta and will expand nearly spherically. Assuming 
that the material inside the cone of a jet-half opening angle 
9j is shocked by the jet and forms a cocoon and the de¬ 
posited energy into the cocoon is E c , the Lorentz factor of 
the cocoon can be estimated as 


1 + 0.05 


Mej(flj) 

10 ~ 5 Mq 


E c 


10 48 erg 


( 1 ) 


where M e j(9j ) is the ejecta mass within 9j. As the jet cross¬ 
ing time is comparable with the duration of a short GRB 
and the jet energy deposited in the cocoon will be compa¬ 
rable to the jet energy, we expect that the cocoon energy 
will be similar to the GRB jet energy. For our fiducial value, 
we take a kinetic energy of 10 48 erg with a single velocity of 
0.3c and the cocoon is sub-relativistic. Note that it becomes 
relativistic at energies higher than 10 49 erg. 


3 THE RADIO SIGNATURE 

The various components of ejecta interact first with each 
other and ultimately with the ISM. This last interaction 
produces a long-lived blast wave. This shock that propa¬ 
gates into the ISM will enhance magnetic fields and accel¬ 
erate electrons that will emit synchrotron radiation. The 
process is similar to GRB afterglows and to radio emission 
from some early supernova remnants. In this section, we ex¬ 
plore the synchrotron radiation from a merger taking into 
account the various components of ejecta. Except for the 
relativistic jet we consider all components of the ejecta as 
spherically symmetric. We discuss the implications of this 
approximation in Sec. m 

We assume that the ISM is homogeneous and charac¬ 
terized by an external density n. The ejecta slows down with 
the deceleration timescale given by 

_ 3 E _ 

47rm p c s nro(ro - 1)/3 q 

where To and /3o are the initial Lorentz factor and the cor¬ 
responding initial velocity of the ejecta, m v is the proton 
mass. The values of fd ec for the different components of 
ejecta are shown in the left panel of Fig. [T] For a mildly or 
sub-relativistic outflow the deceleration timescale character¬ 
izes the observed peak time. For an ultra-relativistic beamed 
jet with a viewing angle 9 a b s > 9j , we have an orphan after¬ 
glow. Namely, we do not see the highly beamed burst and 
early afterglow. But we see the late afterglow when it slows 
down and its less beamed emission includes our line of sight. 
As this happens when V ~ 9 ~ hs . the peak time in the source 
frame is around fdec given by Eq. © with r ~ 0 o bs instead of 
the initial Lorentz factor. Note that for a relativistic outflow 
the observer time is different from the time in source frame 
and it is smaller by a factor of F -2 . Note that for an ob¬ 
served GRB the peak in the radio arises when the observed 
frequency equals the typical synchrotron frequency. 

The deceleration timescale (see the left panel of Fig. [l]) 
suggests three types of the radio flares. First, the ultra- 
relativistic jet produces at early times the radio afterglow, 
that can be seen only by observers along the jet axis or 
close to it. Second, the mildly relativistic components, in¬ 
cluding the cocoon, the shock-breakout material, and the jet 
for an observer away from its axis produce radio flares with 




a timescale of a few dozen days. Finally, the sub-relativistic 
dynamical ejecta produces a late-time radio flare with a 
timescale of a few years. 


3.1 Ultra-relativistic beamed jet 

The Blandford-Mckee self-similar solution describes the jet 
dynamics, in the relativistic regime after the energy of the 
ISM swept up by the jet becomes comparable to the energy 
of the jet itself. Once the Lorentz factor of the jet decreases 
to F ~ fU 1 , the jet expands laterally and approaches a quasi- 
spherical shape. To describe the evolution during this side¬ 
way expansion phase we ad opt a semi-analyt i c form ula for a 
homogeneous jet given by iGranot fc Piranl d2012H 2 l which 
shows a good agreement wit h the results of a numerical 
simulation by iDe Colle et ali (|2012l 'l. The observed signal 
depends strongly on the viewing angle and we consider five 
different viewing angles 9 a b s = (0°, 30°, 45°, 60°, 90°). 

To calculate the synchrotron radio emission we assume, 
as common in GR B afterglows a nd in radio supernovae mod¬ 
eling (see e.g., ISari et al.1 11998I L that the shock generates 
magnetic fields and accelerates electrons to a power law 
distribution N( 7 ) oc 7 _p , where 7 is the Lorentz factor 
of an accelerated electron. The value of p is estimated as 
p a 2.1 - 2.5 in late GRB afterglows and afterglows of 
low luminosity GRBs and as p ~ 2.5 - 3 in typical ra¬ 
dio supernovae (IChevalieiill998h . We assume p = 2.5. The 
total energy of the electrons and the magnetic field inten¬ 
sity are characterized by equipartition parameters: t e and 
tB that are the conversion efficiency from the internal en¬ 
ergy of the shock into the energy of the accelerated electrons 
and magnetic fields, respectively. We set these parameters as 
e e = tB = 0.1. These values are consistent with those eval¬ 
uated from late radio afterglows in long GRBs (Frail et al. 
l200d. l2005lh For our purposes the radio emission is always 
below the cooling frequency hence the system has only two 
characteristic frequencies, the synchrotron frequency of the 
“typical” electron and the self absorption frequency. We im- 
plemen t the effect o f the synchrotron-self absorpt i on fo l- 
lowing IGranot et al.1 (|l999bf i ; iRvbicki fc Lightmanl (|l979t l . 
Once we determine the local emissivity we integrate over the 
inte nsity of each line of sight with a n equal arrival time (see 
e.g.. ISarilll998l : IGranot et al.lll99~9aH 8 l . 

3.2 Mildly and sub-relativistic isotropic 
components. 

We briefly discuss the simple analytic estimates of the ra- 
d io signal s for a mildly and sub-relativistic ejecta (see 
IPiran et al.ll2013l for details). The hydrodynamics of a mildly 
and sub-relativistic blast wave with a kinetic energy E and 
an initial velocity /3o expanding into a homogeneous ISM 


2 We adopt the conical model of IGranot fe Piranl J2OI2I1 . The 
difference of the afterglow flux between their different models is 
a factor 2 ~ 3 during the side-expansion phase. 

3 We do not us e th e afterglow library of 
Ivan Eerten fc MacFadvenl (1201 lh which is incorrect below 
the absorption frequency. Above the absorption frequency our 
light c urves are consistent with those of lvan Eerten &; MacFadvenl 

d201lh . 
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with an external density n can be approximately described 
by /3 = j3o until the deceleration time tdec. The dynamics ap¬ 
proaches to the Sedov-Taylor self-similar solution after tdec- 
The synchrotron emission is slow cooling and it is 
strongly suppressed by the self absorption below the self 
absorption frequency: 


1 

\ Va ’ Aec {th~cj 

2 

1 P+4 

( t ^ tdec)? 


I2a{t) = l 

3p —2 

(3) 

1 

l Ua ’ dec 

i P+4 

(t !> tdec)? 



where 


3P+14 2 +P , 2 (P" 1 ) 

- i -—— O I „ I A \ „ I 1 Q ( „ I 1 \ 


^,dec = 1 GHz n^ye|^ e e ^+ 4 0 **™. (4) 


These expressions are valid for v a > u m , where v m is the syn¬ 
chrotron frequency of electrons with the minimum Lorentz 
factor. Here and elsewhere, Q x denotes the value of Q/ 10 x 
in cgs units. 

For v > v a , the peak flux and the peak time can be 
estimated as 


p+1 P+1 _ ^ 5p — 7 

Tpeak,y>i/ a (t dec ) « 0.8 mjy E i9 n~ 2 


(5) 


( Dl 
V 200 Mpc 


(t+s)' 


and 

1 _ 1 _ 5 

+>M*dec) = ^dec ~ 40 day E£ g n 3 /3 0 3 , 


( 6 ) 


where Dl is the luminosity distance to the source and we 
approximate the Lorentz factor as To — 1 ~ /3q in Eq. ©• 
The peak flux and its time depend sensitively on the external 
density, the kinetic energy, and the initial velocity of the 
ejecta in the optically thin regime. 

For v < v a at tdec, the peak flux and peak timescale are 


« 1 i 


Fpeak,p<i/ a « 0.1 mjy E 4 5 9 n 3 e^ _ 1 e^_ 1 
-2 , .. .6 


(7) 


( 


\200 Mpc 


V ( y v 

) \ 150 MHz/ 


and 


+<p»(idec) ~ 200 da y E ls ( 


150 MHz 


\ ^ 

)“■ ( 8 ) 


In the optically thick regime, the peak flux and its timescale 
depend weakly on the external density and they are inde¬ 
pendent of the initial velocity of the ejecta. The dependence 
on the energy is also weaker than in the optically thin case. 

As the velocity distribution is not uniform, we estimate 
the emission from each shell of matter and combined the 
results. For a given distribution of energies as a function of 
velocity, we divide the outflow into shells. An external ISM 
mass, M(R), swept up at a radius R can be associated with 
each shell such that this mass slows down the shells: 

M(R)(cf 3F) 2 = E{> /3F). (9) 

Once we solve the implicit Eq. ©, we determine the ob¬ 
served light curves for each shell. We then combine the con¬ 
tributions of the different shells to obtain the total light 
curve. In the non-relativistic limit, the ejecta dynamics de¬ 
scribed by Eq. @ is c onsiste nt with the self-similar solution 
derived by IChevalierl dl982l l up to a factor of order unity. 
In the relativistic limit and the case of E(^ 0T) = const, it 


agrees with the Blandford-Mackee self-similar solution again 
up to a factor of order unity. 


3.3 Numerical result 

Figure [2] depicts the resulting radio light curves of the differ¬ 
ent components for our fiducial model (see Table [T] for the 
fiducial parameters). We examine three different values of 
the external density n = 0.01 - 1 cm -3 and we present the 
light curves for two observed frequencies 150 MHz (left pan¬ 
els) and 1.4 GHz (right panels) corresponding to the LOFAR 
and the EVLA of radio telescopes. We set the luminosity dis¬ 
tance of the source to be 200 Mpc, which is roughly the sky 
averaged horizon distance of the advanced GW detectors for 
ns 2 mergers. 

The ultra-relativistic jet always arrives first. This on- 
axis emission of the jet, the GRB radio afterglow, is the 
strongest at 1.4 GHz for low external densities n < 0.1 cm' 3 . 
At 150 MHz this GRB afterglow as well as the other rela¬ 
tivistic component, the shock-breakout material, is strongly 
suppressed by self-absorption and it is much weaker. For 
a generic observer, the off-axis orphan afterglow at view¬ 
ing angles of 60° and even at 45° is always subdominant 
compared to the shock-breakout material and the dynam¬ 
ical ejecta. The mildly relativistic component, that arises 
from the shock-breakout material, peaks later at around 20- 
100 days depending on the observed frequency and external 
density. Finally the sub-relativistic dynamical ejecta arises 
at late times (typically 1000 days). It is always the brightest 
at 150 MHz and it is also brightest at 1.4GHz for higher ex¬ 
ternal densities. For our fiducial parameters, that are based 
on a weak GRB, the radio emission from the sub-relativistic 
cocoon is always negligible. 

As mentioned earlier, at early times, synchrotron self¬ 
absorption strongly suppresses the radio flux at 150 MHz. 
As a result, the peak flux is only F v ~ 0.01 mjy for the 
relativistic components such as the shock-breakout material 
and the off-axis GRB jet for all the cases. As expected from 
Eq. 0. in this case the peak flux depends only on the kinetic 
energy among the parameters of ejecta. Indeed, the dynam¬ 
ical ejecta is the brightest as F v ~ 1 mjy for n > 0.1 cm' 3 
and its peak time is ~ 1000 days. For very low densities 
n < 0.01cm -3 , the on-axis GRB afterglow is comparable to 
the dynamical ejecta flare, peaking at about 20 days with 
F v ~ 0.1 mjy. 

At 1.4 GHz, there are the early and late-time radio 
flares. The relativistic components such as the GRB after¬ 
glows and the shock-breakout material contribute the flare 
at early times as expected from Eqs. m and ©. For low 
densities n < 0.1 cm -3 , the GRB afterglows within a viewing 
angle 0 o b s ~ 30° is the brightest at 1.4 GHz as ~ 0.1 
0.5 mjy. Note that the off-axis GRB afterglows are very 
faint for large viewing angle 0 O b s > 60° compared with the 
shock-breakout material and the dynamical ejecta. At this 
stage the originally beamed jet has already slowed down 
and its emission is already quasi spherical because of its low 
Lorentz factor (this is independent of the question how much 
did the jet physically expand sideways). As the jet energy is 
smaller than those of the other components, its radiation is 
weaker. The dynamical component dominates at late times 
t > 100 days and has a relatively flat light curve. 

As mentioned earlier, our fiducial GRB was a typical, 
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150MHz, n=1.0cm' 3 


1.4GHz, n=1.0cm" 3 




150MHz, n=0.1cm" 3 


1.4GHz, n=0.1cm' 3 




150MHz, n=0.01cm’ 3 


1.4GHz, n=0.01cm" 3 




Figure 2. Radio light curves from the different ejecta components for the fiducial model. Here the energy of the ultra-relativistic jet 
and the cocoon is set to be 10 48 erg and the velocity of the cocoon is 0.3c. The external density and the luminosity distance to the source 
are set to be 0.01 1 cm -3 and 200 Mpc. The left and right panels show the radio light curves at 150 MHz and 1.4 GHz, respectively. 


low luminosity one. Figure [3] shows the radio light curves for 
the case of a strong GRB with a jet energy of 10 49 erg (cor¬ 
responding to an isotropic equivalent energy of ~ 10 51 erg). 
The energy of the cocoon is also larger as this should be 
comparable to the jet energy. We use E c = 10 49 erg and a 
corresponding Lorentz factor of T = 1.5 is obtained from 
Eq. ©. Now, for this GRB, the on-axis GRB afterglow is 
the brightest at all densities at 1.4 GHz and for very low 
densities at 150 MHz. The cocoon is much brighter than 
that of the weak GRB and its peak flux at 1.4 GHz is com¬ 


parable to that of the dynamical ejecta and to the off-axis 
orphan afterglow for 0 O b s = 45°. 

Unfortunately, there are numerous uncertainties in our 
estimated light curves. One of the strongest sources of the 
uncertainties arises from lack of precise estimates of the mass 
and energy of the different components. To demonstrate the 
possible variability of the light curves with these unknown 
parameters, we present in Fig. [I] the dependence of the ra¬ 
dio flares on the kinetic energy of the different components. 
Here we show the light curves of the fiducial model (thick 
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150MHz, n=1.0cm' 3 


1.4GHz, n=1.0cm" 3 




150MHz, n=0.1cm 


-3 


1.4GHz, n=0.1cm' c 




150MHz, n=0.01cm' c 


1.4GHz, n=0.01cm' 3 




kinetic energy of the ultra-relativistic jet and the cocoon is set to be 


Figure 3. Same as Fig. [2] for the case of a strong GRB. The 
10 49 erg. The Lorentz factor of the cocoon is 1.5 (see Eq. [lj. 

curves) and those with kinetic energies larger and smaller by 
a factor of 3 (thin curves) than those of the fiducial model. 
Above (below) the self-absorption frequency, the amplitude 
of the light curves scales as oc E ( E 4,/5 ) and the timescale of 
the light curves behaves as oc E 1 ^ 3 ( E 5 ' n ) as expected (see 
Sec. 13.211 . 

While the overall peak luminosity and peak flux depend 
mostly on the global properties of the outflow component, 
the details depend also on the spatial distribution and on 
the velocity distribution. For instance, the light curve of an 
off-axis jet rises steeply because of the collimation of the jet 


and relativistic beaming effect. The detailed shapes of this 
rise will depend on the angular structure of the jet. The light 
curves of the spherical components rise slowly compared to 
that of off-axis afterglows. The slope of the light curve de¬ 
pends on the velocity distribution E (^ T/3) oc (r/3) -Q . The 
rise of light curves will be shallower for lower values of a. 

Examination of the numerical simulations reveals that 
the mildly and sub-relativistic components that we have ex¬ 
amined do not satisfy indeed the spherical symme t ry as - 
sumption that we have made here. iMargalit fe Piranl (120151 ) 
have estimated the effects of a-sphericity on the emission, 
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150MHz, n=1.0cm' 3 


1.4GHz, n=1.0cm" 3 




150MHz, n=0.1cm" 3 


1.4GHz, n=0.1cm' 3 




150MHz, n=0.01cm’ 3 


1.4GHz, n=0.01cm" 3 



t [day] 



Figure 4. Uncertainties in the radio light curves. The thick lines denote the light curves of the fiducial case (see Fig.(2j and the thin 
curves denote the light curves with a kinetic energy larger and smaller by a factor of 3 than those of the fiducial models. 


focusing on the dynamical ejecta. They found that, for a 
given total mass, energy, and external density, a-sphericity 
typically delays the peak emission and reduces the peak 
flux. This can be understood intuitively as follows. If more 
mass and energy are concentrated in one direction the mat¬ 
ter propagating in that direction will slow down later. This 
longer deceleration time results in a longer and weaker radio 
flare compared with the isotropic one. Note that, however, 
as the outflow is only mildly relativistic, even from a highly 
a-spherical ejecta, the emission will be roughly isotropic and 
viewing angle effects will be small. It is worth noting that the 


effect of a-sphericity is more relevant for black hole neutron 
star mergers, wh ich can result in highly a-spherical mas s 
ejection (see e.g.. lKvutoku et aJ1l2013l : |Foucart et al.ll2013l j. 


4 THE RADIO SIGNATURE OF THE SHORT 
GRB 130603B 

The sho rt GRB 130603B h a d an associated ma cronova can¬ 
didate llBerger et alJ 120131 ; iTanvir et al.l l2013l l. While the 
macronova identification is based only on one observed 
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GRB 130603B, z=0.356 



Figure 5. Radio signatures of the short GRB 130603B and 
light curves at 6.7 GHz (red curves), 1.4 GHz (green curves), 
and 150 MHz (blue curves). The red dotted curve denotes the 
GRB radio afterglow from a jet with E = 8 x 10 48 erg, Qj = 4°, 
p = 2.3, e e = 0.2, cb = 8 x 10 3 , and n = 1.0 cm -3 . The solid 
and dashed curves denote the expected radio light curves from 
a dynamical ejecta with an external density n = 1.0 cm -3 and 
0.5 cm -3 , respectively. For the dynamical ejecta, the kinetic en¬ 
ergy is assumed to be 8 x 10 50 erg and the other microphysics 
parameters are the same as those in Sec. [3] The filled squares and 
the open triangles show the observed dat a points and the upper 
limits at 6.7 GHz obtained with the VLA (iFong et alj|20l3) . The 
blue shaded region shows the expected sensitivity of the EVLA 
at 1.4 GHz. 

data point in the TT-band at about 7 days (in source 
frame) after the burst, this is the first, even though rather 
weak, evidence of the significant mass ejection from a ns 2 . 
Using this data point, one can estimate, from the ob¬ 
served luminosity of the macrono va, the minimal ejecta 
mass as Me\ ~ Q .Q2(eth/0-5) M© (|Hotokezaka et all 120131 ; 
Piran et alj 120141 ), where eth is the conversion efficiency 
from the total energy generated by radioactive decay into 
the thermal energy of the ejecta. The velocity can be 
also estimated as v > 0.1c from the condition that radi¬ 
ation can diffuse out from the ejecta with a mass M e \ > 
0.02M|^ and an opacity o f 10 cm 2 /g dKasen et al.l [201a 
iTanaka fc Hotokezakall2013l ) at about 7 days after the burst. 
Assuming that the ejecta mass and average velocity is M e j ~ 
0.02Mq and v ~ 0.2c, the estimated kinetic energy is about 
10 51 erg. 

Figure [5] shows the expected radio flares from the dy¬ 
namical ejecta at 150MHz, 1.4 GHz, and 6.7 GHz as well 
as the observed early radio afterglow of GRB 130603B at 
6.7 GHz by the VLA (IFong et alj|2014h . Here we also show 
a GRB afterglow light curve which is obtained with pa¬ 
rameters E = 8 x 10 48 erg, 9j = 4°, p = 2.3, e e = 0.2, 
eg = 8 x 1CU 3 , and n = 1.0 cm' 3 . The light curve is consis¬ 
tent with the observed data points and upper limits|3 within 
a factor of 2. It is worth emphasizing that the parameters 
can change by orders of magnitude and still fit the data. For 
instance, the external density lies in the range n « 5 x 10 -3 

4 The light curve that we obtain is also consistent with the ob¬ 
served data in other frequencies except for the late time excess in 
the R-band (the macronova candidate) and in the X-ray band. 


30 cm -3 (IFong et al.l 120141 ). For modeling the radio flare 
from the dynamical ejecta, we assume an external density to 
be 1.0 and 0.5 cm -3 and we use an ejecta mass M = 0.02 M@ 
and a kinetic energy E = 8 x 10 5 ° erg. The predicted dy¬ 
namical ejecta light curves at 6.7 GHz are well below the 
upper limits F„ < 30 /rjy at ~ 80 days. However, later ob¬ 
servations may detect a signal. Specifically the peak flux at 
1.4 GHz can be as high as F v « 20 fi Jy (depending on the 
external density). The expected sensitivity of the EVLA at 
1.4 GHz is also shown in the figure. For the higher range 
of external densities n > 0.5 cm" 3 , the radio flare might 
be detectable with the EVLA. The signal at 150 MHz can 
be F u ~ 30 /rJy, which might be also detectable with the 
LOFAR. A positive detection of a varying radio signal will 
confirm the identification of this event as a ns 2 merger and 
will establish the observed infrared bump as a macronova. 


5 CONCLUSION AND DISCUSSION 

A ns 2 merger ejects a significant amount of mass in several 
different components: a dynamical ejecta, a shock-breakout 
material, a wind from a black hole/neutron star surrounded 
by an accretion disk, and a relativistic jet, producing a GRB. 
As a result of the interaction of the relativistic jet with the 
earlier and slower ejecta along the rotational axis, a cocoon 
is expected to be formed and expands nearly spherically. 
Among the different components of the ejecta, the dynami¬ 
cal ejecta, which is also the most robustly found in numer¬ 
ical simulations, has the largest amount of kinetic energy 
up to E ~ 10 51 erg. This is comparable to the “isotropic 
equivalent” energy of the highly beamed GRB jet. Some¬ 
what surprisingly the beamed GRB jet is among the least 
energetic. Overall we have three important distinct compo¬ 
nents, the dynamical ejecta, which is mostly sub-relativistic, 
the mildly relativistic shock-breakout material and cocoon 
0, and the ultra-relativistic beamed jet. 

We have calculated the expected radio signals pro¬ 
duced via synchrotron emission from electrons accelerated 
in the shocks formed between the different components of 
the ejecta and the ISM. This would be a low frequency 
(as opposed to the X-ray or optical GRB afterglow or the 
optical-infrared macronova) electromagnetic counterpart of 
the GW event. This process is similar to GRB afterglows 
and radio emission of some early supernova remnants. In 
contrast with the high frequency counterparts, this emis¬ 
sion lasts much longer and may even peak a few years after 
the merger. We focused on the expected radio flux at two 
frequencies of 150 MHz and 1.4 GHz. We found that there 
are three types of the radio flares: (i) The ultra-relativistic 
jet produces the earliest bright radio flare with a timescale 
of ~ 10 days for an observer along the jet axis or close to 
it. (ii) A radio flare with a timescale of a few dozen days is 
produced by the mildly relativistic components such as the 
shock-breakout material, the off-axis jet, and the cocoon. 
The latter flare from the cocoon is significant only when 
the cocoon is energetic enough so that it has a relativistic 
velocity, (iii) Finally the sub-relativistic dynamical ejecta 

5 Note that the velocity of the cocoon depends on the energy 
deposited into the cocoon. The cocoon will be relativistic when 
the deposited energy is larger than about 10 49 erg. 
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produces a radio flare with a timescale of a few years. At 
150 MHz, however, the radio flare is strongly suppressed by 
synchrotron self absorption for n > 0.01 cm -3 until a few 
months. Hence the earlier signals are much weaker at low 
frequencies. Depending on the external density and the dis¬ 
tance to the source, these radio flares could be detectable as 
GW counterparts. Such a detection can reveal the nature of 
ns 2 mergers. 

Although the modeling of radio emission from merg¬ 
ers contains uncertain parameters such as the kinetic en¬ 
ergy of the ejecta and the external density, it is worth to 
estimate these values from the nature of the short GRB 
130603B. The detection of a macronova candidate associ¬ 
ated with this event allows us to estimate the ejecta mass 
M e j « 0.02(e t h/0.5) Mg. Assuming the velocity of the ejecta 
is 0.2c, the estimated kinetic energy is about 10 51 erg. The 
detection of the afterglows implies that the ex ternal den¬ 
sity i s in the range of n ~ 5 x 10~ 3 - 30 cnR 3 (iFong et al.i 
120141 ) . The radio afterglow of GRB 130603B, that arose from 
the relativistic jet, decayed quickly and it was below 30 /rjy 
at ~ 4 days with a similar upper limit at 80 days. How¬ 
ever it is still possible to observe the late-long lasting radio 
signal arising from the dynamical ejecta. This signal could 
be as high as 20 /rJy at 1.4 GHz depending mostly on the 
external density. For the higher range of external densities 
n > 0.5 cnC 3 , this would be detectable at 1.4 GHz with 
the EVLA. At 150 MHz, the expected flux is about 30 /rjy, 
which depends weakly on the external density and peaks at 
late times. This flux might be detectable with the LOFAR. 
While a detection is uncertain, a positive radio signal will 
confirm the identification of this event as a ns 2 merger and 
will establish the observed infrared bump as a macronova. 
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